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SECTION  I 
INTRODUCTION 

A  simple  empirically  baaed  mathematical  model  for  constant 
amplitude  loading  fatigue  crack  growth  rate  (PCGR)  test  data 
is  very  useful  for  predicting  the  crack  growth  rate  for  a 
particular  material  at  a  condition  where  test  data  are  non¬ 
existent,  In  this  manner  the  necessity  for  generating  data 
at  a  particular  unexamined  test  condition  is  circumvented. 

There  are  several  models  already  in  existence  that  vary  in 
their  degree  of  complexity  and  their  degree  of  success  in  pre¬ 
dicting  test  data  results.  Reference  1  discusses  a  simple 
empirical  model  for  the  shift  in  the  linear  region  of  room 
temperature  PCGR  data  for  aluminum  alloy  7010-T73651  with  a 
change  in  load  ratio,  R-ratio  (minimum  load/maximum  load). 

The  linear  data  region  of  PCGR  test  data  is  depicted  in 
Figure  1.  The  model  was  based  on  the  Paris  equation; 

da/dn  »  CAk"*  (1) 

where  da/dn  is  the  crack  extension  per  load  cycle,  termed  the 
fatigue  crack  growth  rate,  AK  is  the  stress  intensity  range,  and 
C  and  m  are  material  dependent  constants.  The  Paris  equation 
is  applicable  to  the  linear  data  region  only  (assuming  the  log- 
stress  intensity  range,  log-AK,  is  plotted  versus  the  log- 
crack  growth  rate) ;  the  threshold  and  rapid  growth  rate  regions 
are  not  considered  in  this  expression  (Figure  1) .  The  Reference  1 
model  represents  the  log-Paris  coefficient,  log-C,  as  linearly 
related  to  the  changing  R-ratio  at  room  temperature,  assuming 
a  fixed  exponent.  See  Figure  2  .  The  model  was  successful 
at  accurately  predicting  the  best  fit  straight  line  to  the 
linear  data  region  prior  to  the  generation  of  the  data. 

This  program  expands  the  model  developed  in  Reference  1 
to  account  for  variation  of  the  teat  temperature.  To  accom¬ 
plish  this,  three  issues  will  be  addressed.  At  elevated 
temperatures  does  there  still  exist  a  linear  relationship 
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Figure  2.  Log-Paris  Coefficient  Versus  R-Ratio  for  Aluminuiri  Alloy  7010-T73651 


between  the  R-ratio  and  the  log-Paris  coefficient?  If  so,  is 
there  a  trend  in  the  lines'  slopes  with  a  change  in  test 
temperature?  Lastly,  with  the  R-ratio  held  constant  can  a 
simple  mathematical  relationship  be  defined  to  accommodate  the 
shift  in  the  data's  linear  region  with  a  change  in  test 
temperature? 


SECTION  II 

TEST  PROGRAM  AND  SPECIMENS 


The  test  matei..i.al  was  aluminum  alloy  2024,  half  inch 
(12.7  mm)  thick,  bare,  rolled  plate.  It  was  produced  by  the 
Aluminum  Company  of  America.  The  material  was  provided  in 
the  T351  condition  which  is  a  solution  heat  treatment  followed 
by  cold  workinq  and  natural  aging.  The  results  of  a  chemical 
constituent  analysis  is  presented  as  follows. 

Chemical  Constituent  Composition 

Mg  Mn  Fe  Si  Ti  A1 

1.5  0.58  0.20  0.16  <0.03  Balance 

Tensile  specimens  were  machined  from  the  test  plate  and 
triplicate  tensile  tests  were  performed  at  the  four  test 
temperature  of  interest!  72®,  200®,  300®,  and  400°F  (22°, 

93°,  149®,  and  204°C,  respectively).  The  specimens  were 
machined  in  accord  with  Figure  3-  All  tensile  specimens  were 
fabricated  with  the  loading  direction  parallel  to  the  plate's 
longitudinal  grain  direction.  All  of  these  tests  were  con¬ 
ducted  in  compliance  with  the  applicable  ASTM  test  standard, 

E-8,  "Tensile  Testing  of  Metallic  Materials." 

Two  to  six  constant  amplitude  loading  FCGR  tests  were 
completed  at  each  of  20  different  test  conditions.  The  test 
conditions  were  the  combination  of  five  different  R-ratios: 

0.01,  0.1,  0.3,  0.5,  and  0,6,  and  the  four  different  test  tem¬ 
peratures:  72°,  200°,  300°,  and  400°F  (22°,  93°,  149°,  and 

204 °C,  respectively).  All  of  the  FCGR  tests  were  conducted  in 
accord  with  ASTM  test  procedure  E647-78,  "Constant-Load-Amplitude 
Fatigue  Crack  Growth  Rates  Above  i0~®  m/cycle."  Also,  all  of 
these  FCGR  tests:  (1)  were  conducted  in  a  laboratory  air  en¬ 
vironment,  (2)  used  a  loading  frequency  equal  to  20  Hz,  and 
(3)  used  the  CT  specimen  shown  in  Figure  4  with  L-T  grain 
orientation . 
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Figure  3.  Tensile  Test  Specimen. 
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Figure  4.  CT  Fatigue  Crack  Growth  Test  Specimen. 


Only  the  linear  region  of  the  crack  growth  rate  data 
as  represented  on  a  log-stress  intensity  range  versus 
log-crack  growth  rate  pair  of  axes  was  considered  in 
this  effort;  the  threshold  and  rapid  crack  growth  rate  region 
immediately  preceding  failure  were  not  considered  and  are 
not  presented  herein.  The  maucimum  and  minimum  crack  growth 
rates  that  were  used  to  define  the  linear  region  for  each  of 
the  20  test  conditions  are  listed  in  Table  1;  these  limits 
represent  conservative  subjective  judgements;  generally,  the 
linear  region  extends  beyond  these  limits. 

For  the  first  23  teats  completed,  which  represent  each 
R-ratio  in  combination  with  either  a  72®F  (22®C)  or  200*F 
OS^C)  test  temperature,  the  Paris  exponent  and  coefficient 
were  allowed  to  freely  vary  when  calculating  the  best  fitting 
linear  equation  to  these  individual  specimen's  data  sets. 

The  average  value  exponent  of  these  23  individual  specimen  data 
sets  was  3.36  with  the  maximum  value  of  3.50  and  a  ndnimum  value 
of  3.27  or  a  range  of  plus  or  minus  4  percent.  Subsequently 
in  calculating  the  best  fitting  equation  to  the  20  multi¬ 
specimen  data  sets  the  Paris  exponent  was  fixed  equal  to 
^^=3.36,  and  only  the  Paris  coefficient,  C,  was  allowed  to 
freely  vary. 

Following  the  calculation  of  the  best  fitting  quation 
to  each  of  the  20  data  seta  in  accord  with  the  above  described 
procedure  a  mathematical  model  of  the  shift  in  the  Paris 
coefficient  for  a  change  in  R-ratio  and/or  test  temperature 
was  formulated.  The  formulated  mathematical  model  was  used 
to  predict  the  best  fitting  equation  to  a  test  case  set  of 
data  prior  to  the  generation  of  the  test  case  data.  The  test 
case  was  arbitrarily  selected  to  be  250®F  (121®C) 
at  a  load  ratio  of  0.35. 
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SRCTION  III 

RESULTS  AND  DISCUSSION 

Tensile  test  results  are  presented  in  Table  2.  The 
material  is  a  moderate  strength  and  ductile  aluminum  alloy/ 
heat  treatment.  For  the  temperature  rise  from  room  temperature 
to  400*F  (204®C)  the  average  ultimate  strength  decreases 
27.0  percent,  whereas,  the  average  yield  strength  only  decreases 
10.5  percent.  For  the  same  temperature  rise  there  is  little 
change  in  the  percent  elongation  at  failure,  however,  there  is 
a  large  increase  in  the  percent  reduction  of  area. 

The  linear  region  FCGR  test  results  are  presented  in 
Appendix  A  in  Figures  A.l  through  A. 20.  Generally,  the  linear 
region  .shifts  down  and  to  the  left  with  an  increa.se  in  R-ratio 
and  for  this  material  changes  very  little  with  an  increase  in 
temperature.  The  best  fitting  equation  that  was  calculated 
for  each  data  set  with  the  exponent  fixed  equal  to  m=a3.36  is 
also  listed  on  each  of  the  20  figures.  From  this  point  on  in 
the  discussion  of  modeling  the  FCGR  data  with  the  Paris  equation, 
the  stress  intensity  range,  AK,  is  in  KSl/In,  the  crack  growth 
rate,  da/dn,  is  in  in. /cycle  and  the  temperature  is  in  degrees 
Fahrenheit  (“F). 

Table  3  lists  the  logarithm  of  the  Paris  coefficient, 
log-C,  for  all  20  test  conditions  of  interest.  Figure  5  presents 
a  plot  of  the  loading  ratio  versus  the  log-Paris  coefficient. 

Here  if  the  points  for  a  load  ratio  equal  to  0.01  are  excluded 
the  load  ratio  versus  log-Paris  coefficient  can  fairly  well  be 
represented  as  a  straight  line,  as  was  done  in  Reference  1.  This 
is  true  not  only  of  the  room  temperature  tests  but  is  equally 
applicable  to  the  elevated  temperature  data.  The  lowest  load 
ratio  Paris  coefficients  do  not  coordinate  well  with  the 
coefficient  associated  with  larger  R-ratios.  Similarly,  in 
Reference  1  the  Paris  coefficient  of  some  data  for  load  ratios 
less  than  or  equal  to  zero  was  not  linearly  related  to  those 
coefficients  for  the  same  material  generated  at  higher  R-ratios. 
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TABLE  2 


TABLE  3 

LOG-PARIS  COEFFICIENT  FOR  A12024-T351 


TtiSt 


Temperature 
"F  (*C) 

R-  0.01 

0.1 

^.•1 

0.5 

0.6 

72 

(22) 

-0.613 

-8.447 

-8.350 

-8.277 

-8.231 

200 

(93) 

-8.565 

-B.409 

-8.336 

-8.257 

-8.212 

300 

(149) 

-8.500 

-8.402 

-8.323 

-8.257 

-8.178 

400 

(204) 

-8.513 

-8.356 

-8.275 

-8.205 

-8.145 
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In  general,  there  is  very  little  shift  in  the  data  points 
in  going  from  the  72»F  (22"C)  to  the  400»F  (204®C)  data. 

This  is  a  very  desirable  material  characteristic  but  presents 
a  hurdle  in  this  effort  to  characterize  a  material's  response 
to  a  change  in  temperature. 

In  Figure  5  if  lines  were  drawn  through  the  data  points 
for  each  temperature,  excluding  all  of  the  log-Paris  coefficients 
for  an  R-ratio  equal  to  0.01,  the  slopes  of  the  four  lines  (not 
shown)  would  be  0.417,  0.392,  0.424,  and  0,418  for  the  four 
test  temperatures:  72^F  (22'’C)  ,  200®F  /  300®F  (149“C)  , 

and  400*^  (204“C),  respectively.  Since  (1)  there  is  no  trend 
in  the  slopes,  and  (2)  the  slope  for  72“F  (22®C)  and  400'’P  (204®C), 
the  minimum  and  maximum  test  temperatures,  are  practically 
identical,  the  slope  was  assumed  to  be  constant  over  the  tem¬ 
perature  range  and  is  approximately  equal  to  the  average  of 
the  four  values  E  a  0.413.  Assuming  the  slope  of  the  line 
remains  constant  will  accommodate  considerable  simplification 
of  the  mathematical  model  for  the  FCGR  test  data  since 

log-C  =  constant  +  b*R 

is  equally  applicable  for  any  temperature.  Also  the  Paris 
exponent,  m,  is  assumed  to  be  constant,  m»3.36,  over  the  entire 
R-ratio,  temperature,  and  crack  growth  rate  range  (Table  1) 
included  in  this  program,  which  represents  another  convenient 
simplification. 

Figure  6  presents  the  log-Paria  coefficients  listed 
in  Table  3  along  with  the  temperatures.  For  all  five  R-ratios 
there  is  very  little  change  in  crack  growth  data  (Figures  A.  4 
to  A. 20)  or  in  log-C  with  an  increase  in  temperature.  Here 
again  it  can  be  seen  that  all  of  the  log-C  values  associated 
with  a  loading  ratio  equal  to  0.01  plot  disproportionately  low. 

A  linear  relationship  can  quite  adequately  represent  the 
change  in  log-C  with  test  temperature  for  all  five  of  the  R-ratios. 


14 


If  the  best  fit  lines  (not  shown)  were  drawn  for  each  R-ratio, 
the  slope  of  each  line  would  be:  0.00031,  0.00026,  0.00021, 
0.00019,  and  0.00027  for  loading  ratios  equal  to  0.01,  0.1,  0.3, 
0.5,  and  0.6,  respectively  (since  the  plot  for  an  R-ratio  equal 
to  0.01  is  disproportionally  low  these  log-coefficients  data 
points  will  again  be  disregarded) .  The  average  slope  for  the 
four  lines  associated  with  the  R-ratios  larger  than  0.01  is 
d=0. 00023.  This  average  slope  is  represented  as  only  an 
approximate  value  for  the  change  in  log-C  with  a  change  in 
temperature.  Assuming  the  slope  of  the  lines  remains  constant 
represents  another  convenient  simplification  of  the  mathematical 
model  for  the  PCGR  data.  Consequently, 

log-C  *  constant  +  3’T 

is  equally  applicable  over  the  R-ratio  range  from  0.1  to  0.6. 

It  has  been  assumed  that  the  trivariant  data  (R-ratio, 
test  temperature,  and  log»C)  can  be  graphically  represented  by 
a  series  of  parallel  straight  lines.  Reference  2  presents 
a  least  squares  method  for  calculating  the  coefficients  for 
trivariant  linearly  related  data  in  the  general  form 

z  *  bq  +  aj^x  +  a2y  (2) 

for  a  given  set  of  n  data  points  where 

_ Nl  -  N2 _ 

[n  Ixi-();xi)  2]  [n  Zyi--(  J:yJ)  ]-tn  ExiYi- (^xi)  (Jly;^)  ]  ^  (3) 

where 

Nl  =  [n  Exi'- (Ix^)  [n  lyiZi-(Eyi)  (Lzi)  )  (4) 

N2  »  [n  l-x^yi- (Exi)  (iiyi)  ]  [n  (Ixi)  (Ezi)  ]  (5) 
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■1 


tn  Sx^zi- (Exi)  ( ]  -  a2[n  S^xiyi- (Ixi)  (lyi)  ] 
“  - TT - = - 


16) 


n  EXj^  -  (Ex^) 


and  a  I 


Ezi  -  a2Eyi,  -  a^Exi 


IQ  ...  - 

for  i  =  1 , 2 , 3 .  .  .n 


n) 


For  this  particular  application 

X  «  R-ratio  =  R 
y  a  Test  temperature  (®F)  =  T 
and  z  =  log-C, 

Since  the  data  for  an  R-ratio  equal  to  0.01  appears 
to  be  disproportionately  low,  only  the  16  data  points  (n=16)  for 
R-ratios  greater  than  or  equal  to  0.1  were  used  as  input  to  the 
above  equations.  The  results  of  the  calculations  were 


aQ  «  -8.503 
ai  =  0.412 
and  32  “  0.00023 

or  log  C  «=  -8.503  +  0.412R  +  0.00023T  (8) 


The  two  coefficients  calculated  in  this  manner  are 
practically  identical  to  the  average  values  for  the  lines' 
slopes  (Figures  3  and  4)  presented  as  approximations  earlier. 
Taking  the  antilogarithm  of  equation  (8)  yields  a  general 
expression  for  the  Paris  coefficient,  C, 

C  =  j^g(-8. 503+0. 412R+0.00023T) 

Therefore,  by  substituting  equation  (9)  and  the  average 
exponent,  m=3.36,  into  equation  (1)  the  final  mathematical 
model  for  the  test  material  is; 

da  ,  ,,(-8.503+0.412R+0.00023T)  .„3.36 

=  lU  AK  iiu; 

for  72®  <  T  <  400'=F  and  0.1  <  R  <  0.6. 


I 


17 


A  test  case  to  verify  the  general  Paris  expression  (10)  was 
arbitrarily  selected  at  an  R-ratio  equal  to  0.35  and  a  test 
temperature  of  250®F  (12i'’C)  .  Putting  the  values  of  these  two 
parameters  into  equation  (10)  a  predictive  equation  for  the 
data  is  obtained; 

M  -  in("8-!303+0. 412*0. 35+0. 00023*250)  ,^3.36 
dn  ~ 

^  =  j^q(-8.  302)^j^3.  36 

=  4.99  •  10"^AK^'^®  for  R=0.35  and  T=250‘>F  (121®C) 

Two  speci.Tiens  were  tested  for  the  test  case.  Test  results 
are  presented  in  Figure  7  along  with  the  best  fitting  line 
calculated  with  the  exponent  fixed  equal  to  3.36.  (The  best 
fitting  equation  calculated  for  this  same  data  set  with  both 
the  exponent  and  coefficient  free  to  vary  has  the  same  exponent, 

.Q 

3.36,  and  a  coefficient  equal  to  4.76*10  rather  than 
-9 

4.75*10  .)  Also  presented  in  Figure  7  is  the  predictive 

equation  (11).  The  two  lines  overlap.  They  were  started  and 
ended  at  different  stress  intensity  ranges  to  facilitate  visual 
detection  of  two  different  lines  as  opposed  to  one  broad  line. 
Agreement  between  the  predictive  and  best  fitting  equation  for 
the  test  case  data  set  is  excellent. 
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Figure  7.  R=0.35,  250®F  (121»C)  FCGR  Test  Results. 
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SECTION  IV 
CONCLUSIONS 

The  following  conclusions  are  applicable  for  the  test 
conditions  used  throughout  this  report;  i.e.i  72“F  (22®C)  < 
Temperature  <  400"?  (204*C) ,  0.1  £  R-ratio  £  0.6,  and  a  loading 
frequency  equal  to  20  Hz. 

1.  The  log-Paris  coefficient  can  be  modeled  as  a  linear 
relationship  of  the  R-ratio  and  the  test  temperature 
assuming  all  other  test  parameters  remain  constant. 

2.  In  a  log-stress  intensity  range  versus  the  log-crack 
growth  rate  plot  the  linear  region  shifts  down  and  to 
the  left  with  increasing  R-ratio. 

3.  For  a  constant  loading  frequency  and  load  ratio 
there  is  very  little  acceleration  in  the  crack  growth 
rate  in  aluminum  alloy  2024-T351  with  an  increase  in 
test  temperature. 

4.  The  crack  growth  model  derived  herein  netted  a  good 
fitting  predictive  equation  to  the  linear  region  of  a 
test  case  data  set. 
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Figure  A. 2.  R=0.1,  72*P  (22*C)  FCGR  Test  Results. 
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Figure  A. 3.  R=0.3,  72"?  (22'*C)  FCGR  Test  Results. 
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Figure  A. 4.  R-0.5,  72"F  (22*0  PCGR  Test  Results. 
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Figure  A. 7.  R»0.1,  200®F  OS^C)  FCGR  Test  Results. 
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Figure  A, 8.  R-0.3,  200®F  (93®C)  FCGR  Test  Results. 
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Figure  A.ll.  R»0.01,  300^  (149'’C)  FCGR  Test  Results. 
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Figure  A.  12.  R=«0.1,  300'’F  (lAQ^C)  FCGR  Test  Results. 
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Figure  A. 11.  R*0.01,  300"F  (149'’C)  FCGR  Teat  Results. 
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Figure  A, 12.  R-0.1,  300"F  (149®C)  FCGR  Test  Results. 
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Figure  A. 13.  R«0.3,  300"?  (149®C)  FCGR  Test  Results. 
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Figure  A.  16.  R=«0.01,  ^OO^F  (204®C)  PCGR  Test  Results. 
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Figure  A. 18.  R=«0.3,  400®F  (204‘’C)  FCGR  Test  Results. 
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Figure  A. 19.  R-0.5,  400*^  (204'’C)  PCGR  Test  Results. 
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